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ABSTRACT 

We present radiative transfer models for submillimeter galaxies with spectroscopic redshifts and mid-infrared spec- 
troscopy from Spitzer/IRS and analyze available Spitzer/MIPS 24, 70 and 160/im data. We use two types of starburst 
models, a cirrus model and a model for the emission of an AGN torus in order to investigate the nature of these objects. 
We find that for three of the objects (25%) cirrus emission alone can account for the mid-infrared spectrum and the 
MIPS and submillimeter data. For the remaining objects we find that we need a combination of starburst and cirrus in 
order to fit simultaneously the multi-wavelength data. We find that the typical submillimeter galaxy has comparable 
luminosity in the starburst (median L — 10^^'^ L©) and cirrus (median L = lO^'^'* L©) components. This could arise 
if the galaxies have been forming stars continuously for the last 250Myr with the star formation occurring in the last 
5Myr being shrouded by high-optical-depth molecular cloud dust, whereas the rest of the starlight is attenuated by 
diffuse dust or cirrus with an Ay of about Imag. 
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1. Introduction 

The discovery of submillimeter galaxies (SMG) with the 
SCUBA instrument mounted on the JCMT about a decade 
ago (Small et al. 1997, Hughes et al. 1998, Barger et al. 
1998), and the realization that most of them are at high red- 
shift (Chapman et al. 2005, Dannerbauer et al. 2004), was 
of particular significance for studies of galaxy formation and 
evolution. The implied high bolometric luminosities of the 

' SMGs prompted suggestions that they are distant analogs 
of local ultra-luminous infrared galaxies (ULIRGs) which 
emit most of their energy in the far-infrared part of the 

' spectrum. 

Over the last decade significant progress has been made 
■ in the understanding of SMGs. Ivison et al. (2002) through 
\ deep radio mapping of the areas covered by the 8mJy 
. SCUBA survey (Scott et al. 2002) were able to identify 
' the radio and optical counterparts of a significant frac- 
tion of the detected SCUBA sources. This allowed the de- 
termination of photometric redshifts for the sources using 
the far- infrared/radio correlation and other methods and 
the first estimates of their luminosities and star formation 
rates which are found to be of the order of 1000 M^lyr. 
Chapman et al. (2005) were able to obtain spectroscopic 
redshifts for the radio identified sources and determine a 
median redshift of 2.4. Alexander et al. (2005) studied the 
X-ray properties of SMGs and concluded that the majority 
host an AGN which, however, is not luminous enough to 
dominate the bolometric luminosity. 

The diagnostic power of mid-infrared spectroscopy was 
first recognized by the pioneering ground-based work of 
Roche et al. (1991). These studies showed that starburst 
galaxies display emission features that are attributed to 
PAH molecules whereas these features were absent from the 
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spectra of active galactic nuclei. ISO studies (e.g. Genzel 
et al. 1998) used mid-infrared spectroscopy to show that 
ULIRGs are primarily powered by star formation. More 
recently Spitzer studies extended this work to large sam- 
ples of infrared galaxies (e.g. Hao et al. 2007, Spoon et al. 
2007). The unprecedented sensitivity of the IRS on-board 
Spitzer made mid-infrared spectroscopy of SMGs possi- 
ble (Menendez-Delmestre et al. 2007, Valiante et al. 2007, 
Pope et al. 2008). The spectra show that SMGs display 
band emission of polycyclic aromatic hydrocarbon (PAH) 
molecules and have similar mid-infrared spectra as the local 
starburst galaxy M82. 

In the local Universe, the infrared emission of nor- 
mal star-forming galaxies can be understood in terms of 
two components (Rowan- Robinson & Crawford 1989): star- 
burst emission which is associated with optically thick gi- 
ant molecular clouds illuminated by recently formed stars 
and cirrus emission which is associated with diffuse and 
cold dust (r < 30K) illuminated by the interstellar radia- 
tion field. At a time when spectroscopic redshifts for SMGs 
were scarce, Efstathiou & Rowan-Robinson (2003; hereafter 
ERR03) showed that good fits to the spectral energy distri- 
butions (SED) of SMGs with radio or millimeter detections 
could be obtained with pure cirrus models or a combination 
of cirrus and starburst with the submillimeter dominated 
by cirrus. The high z cirrus models of ERR03 assumed an 
optical depth and intensity which is a factor of 2-3 higher 
than the values found in local cirrus galaxies and the inter- 
stellar radiation field was assumed to be that of a 250Myr 
old burst of star formation. 

Over the last two decades we have developed 
fairly sophisticated radiative transfer models for AGN 
torus (Efstathiou & Rowan-Robinson 1990, 1995), star- 
burst (Rowan-Robinson & Efstathiou 1993, Kriigel & 
Siebenmorgen 1994, Efstathiou et al. 2000, Siebenmorgen 
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& Kriigel 2007) and cirrus (Siebenmorgen & Kriigel 1992, 
Efstathiou & Rowan- Robinson 2003). The starburst and 
cirrus models take into account the effect of small grains 
and PAH molecules. Starburst models have also been devel- 
oped by Silva et al. (1998), Takagi et al. (2003) and Dopita 
et al. (2005). Our cirrus model does not take into account 
the distribution of dust and stars as the GRASIL (Silva 
et al. 1998) or the Piovan et al. (2006) models do but it 
has been shown by ERR03 to give spectra that are in good 
agreement with the spectral energy distributions of local 
cirrus galaxies. As we show in section 3, the model is also 
in very good agreement with observations of cirrus in our 
own galaxy. Other work on radiative transfer modeling of 
galaxies has been presented by Bianchi et al (1996), Xilouris 
et al. (1999) and Popescu et al (2000). Other work on ra- 
diative transfer modeling of the torus in AGN has been pre- 
sented by Pier & KroHk (1992), Granato & Danese (1994), 
Nenkova et al. (2002, 2008), DuUemond & van Bemmel 

(2005) , Honig et al. (2006) and Schartmann et al. (2008). 
In this paper we use radiative transfer models of star- 
burst, cirrus and AGN torus emission to constrain the 
properties of SMGs with spectroscopic redshifts and mid- 
infrared spectroscopy and far-infrared photometry (Sect. 2) 
from Spitzer. For comparison two different starburst mod- 
els are applied: an evolutionary model (Efstathiou et al. 
2000) that incorporates a stellar population synthesis model 
and therefore provides information about the stellar pop- 
ulation that powers the starburst (Sect. 4) and a hot spot 
model (Siebenmorgen & Kriigel 2007) that is sensitive to 
the geometry of the dust and stars and therefore pro- 
vides information about the size of the starburst region 
(Sect. 5). A flat Universe is assumed with A = 0.73 and 
Ho = 71km/s/Mpc. 

2. Observations 

Our sample is determined solely by the requirement that 
the galaxies have been detected in the submillimeter and 
have mid-infrared spectroscopy from the infrared spec- 
trograph (IRS, Houck et al. 2004) of the Spitzer Space 
Telescope (SST, Werner et al. 2004). We model all the 
galaxies in the PAH sample of Menedez-Delmestre et al. 

(2006) , except SMM J221733+001120 which is at low red- 
shift and which is clearly dominated by cirrus emission, and 
the objects in the sample of VaHante et al. (2007). The data 
in the sample of Menedez-Delmestre et al. cover the silicate 
feature and therefore provide information about the extinc- 
tion in the galaxies but the objects in the sample of Valiante 
et al., being at higher redshift, do not cover completely the 
silicate band. In total there are 12 objects spanning a red- 
shift range from 1.2 to 3.4 (Table 1). Besides IRS, ISOCAM 
(Webb et al. 2003) and submillimeter photometry (Small et 
al., 2002, Scott et al., 2006, Ivison et al. 2005, Zemcov et 
al. 2007) we complement the SED by retrieving all publicly 
available data of the Spitzer far infrared imager (MIPS, 
Rieke et al. 2004) of the targets in the three photometric 
channels centered at 24, 70 and 160/im. 

MIPS raw data are processed by the Spitzer pipeline 
(version S16.1; Gordon et al. 2005) to a flux calibrated mo- 
saic image. The various mosaic images of a particular target 
and channel and from the different programs and observers 
are combined to a final image using SWARFQ. For the 70 

^ available at: |http: //terapix. iap.fr/sof t/swarp| 



Table 1. Sample of SMGs with redshift and MIPS pho- 
tometry or 3cr upper limits (in mjy). 



Name 


z 


24/im 


70pim 


160/im 


SMM J163658.78+405728.1 


1.2 


0.45 


<4.9 


<30 


SMM J030227.73+000653.5 


1.4 


0.23 


<8.8 




MM J163639-I-4056 


1.5 


0.23 


<5.5 


<21 


MM J163650-I-4057 


2.4 


1.02 


<5.5 


<22 


SMM J09429+4659 


2.4 


0.19 


<6.8 


<23 


MM J163706+4053 


2.4 


0.41 


<5.1 


<17 


MM J105155-H5723 


2.7 


0.18 


<2.2 


<16 


SMM J105207.56+571904.7 


2.7 


0.24 


<2.0 


<13 


SMM J00266-H708 


2.7 


0.32 






MM J154127-I-6616 


2.8 


0.24 






SMM J02399-0136 


2.8 


1.24 


15t 


16 


SMM J09431-^-4700 


3.4 


1.04 


<7.2 


<22 



Tentative detection: brightest source pixel appears ~ 8'' off from NED 



position. 



and 160/im bands we use the filtered mosaic images of the 
pipeline (see MIPS Data Handbook). The final image has 
higher redundancy and signal-to-noise than the one ob- 
tained from data of a particular observing sequence. In the 
24/im band all sources are detected and final images are 
shown in Fig. [TJ In this band SMGs appear well separated 
from other sources. In the other channels, at 70 and 160/im, 
with the exception of SMM J02399-0136, SMGs remain un- 
detected. The flux is derived using an aperture centered on 
the first Airy ring and a 2 pixel wide background annulus 
outside the second Airy ring; colour and PSF correction 
factors are applied. The photometric error is better than 
10%. MIPS photometry of the SMGs is given in Table [T] 
and agree with IRS. We verified our procedure on the cal- 
ibration standard star HD106252. For this star aU MIPS 
24/im data are pipeline processed and mosaic images coad- 
ded with SWARP to a final image. On this image we mea- 
sure a flux which is consistent to within 2% with the flux 
measured by Engelbracht et al. (2007) on the same star. 

3. Cirrus and starburst-only models 

We first compare the data of the galaxies in our sample 
with pure cirrus and pure starburst models. For the cirrus 
emission we follow an approach similar to that of ERR03 
except that we do not attempt to link the optical-UV emis- 
sion with the infrared in a self-consistent manner. We first 
determine the near-infrared to UV spectrum by assuming 
an age for the galaxy and a star formation history. We use 
the stellar population synthesis model of Bruzual & Chariot 
(1993) with a Salpcter IMF and stellar masses in the range 
0.1-125 Mq. For this study we assume, as in ERR03, an 
age of 250Myrs and a star formation rate that is constant 
with time. For the last 5Myrs we assume that the stars are 
embedded in the molecular clouds in which they formed 
so they do not contribute to the starlight that is illumi- 
nating the cirrus dust. We then use the parameter to 
scale the spectrum emitted by the stars; "0 is defined to 
be the ratio of the bolometric intensity of starlight to the 
bolometric intensity of the interstellar radiation field in the 
solar neighborhood (Mathis et al. 1983). In Fig. [2] we com- 
pare the spectrum emitted by dust that is illuminated by 
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Fig. 1. MIPS 24/xm images of tlie SMGs listed in Table [T] Images are background subtracted, centered on target and 
displayed in a field-of-view of 100" x 100". Gray scale index in MJy/sr as indicated and contours at 50%, 33%, 25% and 
20% level from the brightest pixel in the image. North is up and East is left. 
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Fig. 2. Comparison of cirrus models with DIRBE observations for High Galactic Latitudes (Arendt et al. 1998). The 
black solid line shows the emission from dust illuminated by the Mathis et al. (1983) interstellar radiation field and a 
•0 of 1. The dust model is described in Sect. 3. The magenta dashed line shows the emission from a dust mixture more 
typical for a starburst where the PAH abundance is reduced by a factor of 5. 
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Fig. 3. Comparison of the spectra of SMGs with pure cirrus (blue dashed) and pure starburst (purple dotted) models. 
The cirrus model assumes that the stellar population that illuminates the dust is 250Myrs old and has a bolometric 
intensity which is '0 = 5 times that of the interstellar radiation field in the solar neighborhood. The starburst model 
assumes a nuclear radius of 3kpc, a visual extinction of 18mag and a luminosity of L = 10^^'^Lobut scaled to fit the 
mid-infrared photometry. Data from this work and Small et al. (2002), Scott et al. (2002), Webb et al. (2003), Ivison et 
al. (2005), Greve et al. (2005), Scott et al. (2006), Menendez-Delmestre et al. (2006), Kovacs et al. (2006), Vahante et 
al. (2007) and Zemcov et al. (2007). For each galaxy two panels are shown: one covering the 3 — 600/im and the other 
the 6 — 12/im rest-frame wavelength range; the vertical scale of the 6 — 12/im panel is linear between — l.SmJy. 
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the Mathis et al. interstellar radiation field for the solar 
neighborhood {ip — 1) with DIRBE data for High Galactic 
Latitudes (Arendt et al. 1998). To determine the spec- 
trum emitted by the dust we use the latest dust model by 
Siebenmorgen et al. (2001). It consists of silicate and amor- 
phous carbon grains with a size distribution: 71(a) oc a^^-^, 
300A < a < 2400A, a population of small graphite grains 
with n(a) oc a''^, lOA < a < 8OA and PAHs. There are 
small PAHs with Nc = 50 and Nh ^ 12 and large PAHs 
with Nc = 300 and Nh = 48, where Nc is the number 
of C atoms and Nh the number of H atoms of a PAH. By 
mass 63% of the dust is in silicates and 37% in carbon of 
which 80% is amorphous, 10% graphitic, 5% in small and 
5% in large PAHs. The element abundance, with respect 
to H, of C and Si in grains is 200ppm and 31ppm, respec- 
tively. The gas-to-dust mass ratio, which allows conver- 
sion from dust into gas mass, is 125. Such a dust mixture 
produces a reddening in rough agreement with the stan- 
dard interstellar extinction curve for Ry = 3.1 and can be 
considered to be typical for interstellar dust in a quiescent 
galaxy like our own. However, in the more active environ- 
ment of a starburst galaxy some of the small grains and 
PAHs are destroyed by the harder photon and stronger ra- 
diation field (e.g. Omont 1986, Leach et al. 1989a, b, Voit 
1991, Voit 1992, Rapacioh et al. 2006, Siebenmorgen & 
Kriigel 2009), thus reducing their abundance. Vega et al. 
(2005) even suggest that all PAHs are destroyed in the 
molecular clouds that constitute the starburst and that the 
PAH emission of infrared galaxies arises from cirrus dust. 
The predicted spectrum of the dust emission that is illu- 
minated by the Mathis et al. interstellar radiation field but 
with reduced small grain abundances is shown in Fig [21 We 
consider that for a starburst 5% of carbon in the grains is 
in small graphites and 2% of C is in PAHs. For the rest of 
the paper we will assume these parameters of the dust for 
our cirrus and starburst model. 

In Fig. [3] we compare our cirrus model with the data 
of the galaxies in our sample. Clearly we cannot constrain 
the value of tp with the presently available far-infrared and 
submillimeter data. However, in order to be consistent with 
the finding of Coppin et al. (2008) that most SMGs emit 
at a temperature of 28K we use ?/; = 5. It is important to 
note that the estimate of Coppin et al. is based on fitting a 
modified blackbody to the 350/im and 850/xm data whereas 
our cirrus model considers a distribution of grain species 
and sizes each with its own temperature. After normalizing 
the cirrus model to the 850/im data point we find that 
for a couple of the objects (SMM J163658. 78+405728.1 at 
a redshift of 1.2 and MM J163639-f 4056 at a redshift of 
1.5) the cirrus model seems to be adequate for explaining 
the spectral energy distribution. For SMM J10521-I-5719 a 
cirrus model with = 18 is also adequate for explaining 
the SED. For the rest of the objects the model falls short 
of matching the observed mid-infrared spectrum suggesting 
either that the PAH abundance is higher than assumed or 
there is contribution in the mid-infrared from a starburst. 
In Sect. 4 and 5 we will explore the latter possibility. 

In Fig. [3] we also explore the possibility that a starburst 
alone can explain the complete spectral energy distribution. 
To do that we take a single starburst spectrum from the 
SED library of Siebenmorgen & Kriigel (2007). We choose 
a starburst model with a total luminosity of 10^^'^Lq, a nu- 
clear radius of 3kpc and a visual extinction from the edge to 
the center of 18mag. For each galaxy this model spectrum is 



normalized to the mid-infrared observations. It is clear that 
a starburst-only model falls short of matching the submil- 
limeter photometry by at least an order of magnitude. The 
model is described in Sect. 5 but it is important to note here 
that as the starburst radius increases the dust gets colder 
and the spectrum shifts to longer wavelengths. Starburst 
models with a radius of lOkpc come close to matching the 
complete spectral energy distribution but as we discuss in 
Sect. 6 such extended starbursts are excluded by millimetre 
interferometry. 

4. Evolutionary starburst models 

Efstathiou, Rowan-Robinson & Siebenmorgen (2000; here- 
after ERRS) presented a starburst model that combined 
a simple model for the evolution of giant molecular clouds, 
the stellar population synthesis model of Bruzual & Chariot 
(1993) and detailed radiative transfer that included the ef- 
fect of temperature fluctuating small grains with dust par- 
ticle radius a < lOOA and PAHs to account for the detected 
infrared emission bands (Siebenmorgen & Kriigel 1992). An 
important feature of the ERRS model is that the expansion 
of the HII region leads to the formation of a narrow shell of 
gas and dust. This naturally explains the fact that the near- 
and mid-infrared spectra of starburst galaxies are not dom- 
inated by emission from hot large (a > lOOA) dust grains 
but by the PAHs emission. As in ERRS we assume that 
the initial Ay of the molecular clouds that constitute the 
starburst is 50mag but we use the grain model described 
in Sect. 3. We further assume a constant star formation 
rate and an age of 5Myr. The assumed starburst age is 
very poorly constrained by our modeling but as we discuss 
in section 6 a value of 5Myr can explain the fact that the 
inferred luinosities of the starburst and cirrus components 
are comparable. 

In Fig. m we present fits to the galaxies in our sample 
with a combination of starburst and cirrus. We first normal- 
ize the cirrus model at 850/Ltm and then scale the starburst 
model so that the combination of starburst and cirrus gives 
the best fit to the mid-infrared spectroscopy and is con- 
sistent with the far-infrared photometry. As discussed in 
Sect. 3, for three of the objects we do not need any star- 
burst contribution for explaining the SED. For two of the 
objects we also find evidence for an AGN component which 
we model with the tapered discs of Efstathiou & Rowan- 
Robinson (1995). The more luminous of the two objects 
(SMM J02399-0136), has been classified as a Seyfert 2 by 
Small et al. (2002). SMM J02399-0136 is also the object 
that has been detected at 70 and 160/im. The luminosities 
of the three components and the associated dust masses 
are given in Table 2. In the ERRS model an increase in 
the luminosity translates into an increase in the number of 
molecular clouds that constitute the starburst and therefore 
the dust mass. 



5. Hot spot starburst models 

Siebenmorgen & Kriigel (2007) presented a starburst model 
which is an evolution of an earlier model described by 
Kriigel & Tutukov (1978) and Kriigel & Siebenmorgen 
(1994). The model assumes that the stars are divided in 
two classes: OB stars that are surrounded by dense clouds 
and constitute so-called hot spots and other stars (old bulge 
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Fig. 4. Comparison of the spectral energy distributions of the SMGs with a combination of evolutionary starburst, cirrus 
and AGN torus models. Unless otherwise indicated the cirrus emission is computed with a ^ = 5 scaling of the interstellar 
radiation field. The data and layout of the figure is as described in FigO Cirrus (blue dashed), AGN torus (light blue 
dash-dotted) and total emission spectrum (magenta solid) is indicated. Derived quantities are given in Table [2] 



stars or massive stars) that are dispersed in the diffuse 
medium (see below). The hot spots determine the mid in- 
frared part of the emission spectrum. The outer radius of 
the hot spots is determined by the condition of equal heat- 
ing of the dust by the stars and the ambient interstellar 
radiation field. Both classes of stars are represented in the 
equation of radiative transfer by continuously distributed 
source terms. It is assumed that the number density of the 
hot spots and of the other stars falls off with the radius of 
the starburst as r^^-^. 

In addition to the dust in the hot spots the model as- 
sumes that the volume of the starburst is filled by dust 
which is uniformly distributed and gives rise to a total ex- 
tinction Ay from the outer radius R of the galactic nucleus 
to its center. In this homogeneous density model the pa- 
rameter Ay is directly related to the dust mass Mgg and 
only one of them is independent. The other model parame- 
ters are the total luminosity Lgg and the starburst radius 
R. The OB stars are assumed to be confined to the central 
350pc whereas the bulge stars fill the whole volume. 

In Fig. [5] we combine starburst models computed with 
the method of Siebenmorgen & Kriigel (2007) with cirrus 
models and compare them with the data of the objects 
in our sample. As in the case of the evolutionary mod- 
els we first normalize the cirrus model at 850fim. Then 
we search in the SED library to find a starburst model 
which, after scaling to the distance of the object but with- 
out further normalization, best fit the mid-infrared spec- 
troscopy. The starburst radius is fixed at 3kpc to be con- 
sistent with the sizes inferred from interferometry for these 
galaxies (Sect. 6). Unless otherwise indicated the value of 
ip is assumed to be 5 as before. 



6. Discussion and Conclusions 

Since the original suggestion of ERR03 and Kaviani et al. 
(2003) that SMGs are colder and more extended than local 
ULIRGs a number of studies have given support to this 
idea. Rowan-Robinson et al. (2004) discovered a number of 
luminous cirrus galaxies in the ELAIS survey. Chapman et 
al. (2004) found that the majority of SMGs in their sample 
are more extended in the radio than local ULIRGs. This 
finding was confirmed by Biggs & Ivison (2008) who found a 
median size of 5kpc for the SMGs in their sample. However, 
Tacconi et al. (2006) showed that the SMGs in their sample 
cannot be more extended than ~ 4kpc. 

The evolutionary starburst model of ERRS does not 
make a prediction about the size of the starburst but a 
lower limit can be obtained by calculating the minimum 
volume needed to contain the molecular clouds that consti- 
tute it. As noted by ERRS a molecular cloud illuminated 
by a lOMyr old instantaneous burst has a luminosity of 
~ 3 X 10® Lq and a radius of 50pc. The minimum radius 
R of the sphere that is needed to contain a starburst of lu- 
minosity L is therefore R = 0.05(L/3 x 10® Lo)!/^ kpc. So 
for a ULIRG with a luminosity of 3 x 10^^ L©, the radius 
is i? ~ Ikpc. A plausible geometry for SMGs is therefore 
a region of 2 — 4kpc radius filled with diffuse dust and in 
which an ensemble of molecular clouds is embedded. 

Kovacs et al. (2006) reported 350/im observations of 
SMGs which suggest that they are colder (mean temper- 
ature ~ 35K) than local ULIRGs. This was confirmed by 
Coppin et al. (2008) who reported 350/xm observations of 
SMGs found in the SHADES survey. In fact Coppin et 
al. find that most SMGs in their sample are actually even 
colder (28K). 
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Fig. 5. Comparison of the spectral energy distributions of the SMGs with a combination of hot spot starburst and cirrus 
models. Unless otherwise indicated all cirrus models assume a "0 = 5 scaling of the interstellar radiation field. Data and 
layout of the figure is as in Fig.O Cirrus (blue dashed) and total emission spectrum (magenta solid) is indicated. Derived 
quantities are given in Table [2l 

Table 2. Derived parameters for each galaxy. 



Name 




Evolutionary models 




Hot spot models 




log LsB 


logic 


log Lagn 


Mc 


log ■'^SB 


log L'q 






[Lo] 


[L©] 


[Lo] 


[IO^Mq] 


[L©] 


[Lq] 


[mag] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SMM J163658. 78+405728.1 




12.2 




1.3 




12.2 




SMM J030227. 73+000653. 5 


11.7 


12.2 


11.3 


1.5 


12.0 


12.3 


7 


MM J163639+4056 




12.3 




1.4 




12.3 




MM J163650+4057 


12.8 


12.4 




2.2 


12.7 


12.6 


9 


SMM J09429+4659 


12.5 


12.3 




1.5 


12.5 


12.4 


36 


MM J163706+4053 


12.4 


12.6 




3.0 


12.4 


12.7 


18 


MM J105155+5723 


12.5 


12.2 




1.3 


12.3 


12.4 


18 


SMM J105207. 56+571904. 7 




12.9 




3.6 




12.9 




SMM J00266+1708 


12.5 


12.8 




4.5 


12.6 


12.9 


36 


MM J154127+6616 


12.4 


12.7 




3.7 


12.3 


12.8 


36 


SMM J02399-0136 


12.7 


12.9 


12.7 


6.0 


13.0 


13.0 


18 


SMM J09431+4700 


13.2 


12.0 




0.9 


13.2 


12.3 


18 



Column (1): starburst luminosity of the evolutionary model. Column (2): cirrus luminosity of the evolutionary model. Column (3): AGN torus 
luminosity. Column (4): dust mass of the cirrus component. Column(5): starburst luminosity of the hot spot model. Column (6): cirrus luminosity of 
the hot spot model. Column (7); visual extinction of the hot spot model measured from the surface to the center of the galaxy. 



Clements et al. (2008) presented spectral energy dis- 
tributions of SMGs detected in the SHADES survey and 
for which there are also supplementary optical and Spitzer 
data. They found that out of the 33 sources in their sam- 
ple 8 can be fitted by a cirrus template whereas most of 
the other sources can be fitted by an Arp220 template. The 
fraction of cirrus dominated galaxies is therefore in good 
agreement with our estimate. The problem with the Arp220 
interpretation for the other objects, however, is that none 



of the objects with mid-infrared spectroscopy from Spitzer 
show the deep silicate absorption features characteristic of 
the Arp220 template. 

Pope et al. (2008) presented mid-infrared spectroscopy 
from Spitzer for a sample of 13 sources in the GOODS 
field. They find that the spectroscopy cannot be fitted with 
an Arp200-type template but instead the average spectrum 
can be fitted by an M82-type one. 
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Farrah et al. (2008) presented mid-infrared spectra of 
high-redshift ULIRGs in a very narrow redshift range of 
1.71 ± 0.15. The ULIRGs in the sample of Farrah et al. 
show mid-infrared spectra similar to local starburst galaxies 
which are two or more orders of magnitude less luminous 
in the infrared. Farrah et al. suggest that one of the most 
likely interpretations of their result is that star formation in 
the ULIRGs of their sample is extended on scales of l-4kpc 
which is similar to the sizes of SMGs observed with radio 
or millimeter interferometry. 

The median dust mass we derive for the cirrus compo- 
nent is 2.2 X 10® M0. For comparison Coppin et al. find a 
median value of 0.9 x 10® Mq but note that most of the 
difference can be attributed to a factor of three difference in 
the assumed absorption coefRcient in the submillimeter. As 
can be seen from Table [2] the luminosities in the starburst 
and cirrus components are comparable. 

The picture we are proposing for SMGs is the follow- 
ing: We assume that stars have been forming continuously 
in the galaxy for the last 250Myr. Our assumption about 
the age of the galaxies is in agreement with the estimate of 
Tacconi et al. (2008). We further assume that the stars that 
formed in the last 5Myr are still embedded in the molecu- 
lar clouds in which they formed. Using Bruzual & Chariot 
(1993) models with a star formation rate of 1000 M0/yr we 
find that the bolometric luminosity of the 5Myr old star- 
burst is 3.2 X 10^^ Lq whereas the luminosity of the rest of 
the stars, which are between 5 - 250Myr old, is 3.5 x 10^^ Lq. 
Assuming that the optical/UV radiation emitted by the 5- 
250Myrs old stellar population is obscured by dust with an 
Ay of Imag, and that all the energy absorbed in the opti- 
cal/UV is re-radiated at a temperature of 28K, this gives a 
cold dust luminosity of 2.5 x 10^^ Lq. So this explains the 
fact that the derived luminosities of the starburst and cir- 
rus components are comparable. The median of the log of 
luminosities of the cirrus component of the 12 galaxies (in 
solar units) is 12.4 for the evolutionary models and 12.6 for 
the hot spot models. The star formation rate for the typical 
galaxy is therefore ~ 1000 — 1600 Mo/yr which is in good 
agreement with the estimates of Ivison et al. (2002). It is 
interesting to note that the timescale we infer for the molec- 
ular cloud phase of star formation is similar to the value 
found by Granato et al. (2000) and Efstathiou & Rowan- 
Robinson (2003) for normal quiescently star-forming galax- 
ies. 

The analysis presented in this paper is clearly limited 
by the lack of data in the rest-frame far- infrared where both 
the starburst and cirrus components peak. Measurements 
with Herschel, whose launch is now imminent, will allow 
us to test the model outlined above. They will allow us in 
particular to put a stronger constraint on the intensity of 
starlight in the case of the cirrus component and possibly 
differentiate between the two starburst models considered 
in this paper. 
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